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Abstract. - 

The solution to the long standing problem of the cohesion of organic chain compounds is 
proposed. We consider the tight-binding dielectric matrix with two electronic bands per chain, 
determine the corresponding hybridized collective modes, and show that three among them are 
considerably softened due to strong dipole-dipole and monopole-dipole interactions. By this 
we explain the unusual low frequency optical activity of TTF-TCNQ, including the observed 
lOmeV anomaly. The softening of the modes also explains the cohesion of the mixed-stack 
lattice, the fractional charge transfer almost independent of the material, and the formation of 
the charged sheets in some compounds. 



Ever since the discovery of the mixed stack organic chain compounds (OCCs), some twenty 
five years ago, two important questions remained open The first is the stability of 

their crystal structure, with a fractional charge transfer between the molecules of n a m 0.6 
electrons per molecule. Furthermore, in some lattices like TTF-TCNQ the equally charged 
(TTF or TCNQ) chains stack side by side, to form charged sheets, instead of alternating in 
both directions, as suggested by the simple Madelung argument. Such a simple alternation 
does occur in other lattices, like HMTTF-TCNQ, which, in spite of the radically different 
chain stacking, have similar n a . Neither the fractional charge transfer nor the formation 
of charged sheets are expected within the standard ionic cohesion scheme which assumes the 
dominance of affinity- ionicity and electrostatic (Madelung) contributions. Besides, calculations 
of cohesion energy within this scheme suggest that almost 2 eV per donor-acceptor molecular 
pair is lacking. This "Madelung defect" [|| clearly signifies that the cohesion at a scale 
of few eVs has to be stabilized by some mechanism(s) beyond the standard ionic scheme. 
Although various suggestions H-H have been put forward, they have failed to explain all 
the aspects of the cohesion, mostly because they could not account for the Madelung defect. 
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On the other hand, Friedel || pointed out that the van der Waals interaction between the 
molecules favors the formation of charged sheets observed in TTF-TCNQ. In this model 
the intramolecular (interband) and the intermolecular (intraband) charge fluctuations are 
separated. Such separation however fails in the low frequency range in which TTF-TCNQ 
shows strong anomalies, observed in numerous microwave/infrared optical measurements 

m (ng.0). 

The origin of these anomalies is the second puzzle of OCCs. E. g. in TTF-TCNQ the 
frequency of the anomaly at 10 meV jl0|, [ll| is comparable to the value of the (pseudo)gap 
produced by the 2fcj? lattice instability at low temperatures (T < 50K). However the anomalies 
in fig.0 can not be related to this instability, since they persist at high temperatures, where 
the 2kp lattice fluctuations are practically absent. 

In this Letter we present a theory which for the first time treats interband (van der Waals) 
and intraband correlations on equal footing, and explains simultaneously the crystal stability 
and the optical anomalies of OCCs. We investigate a simple model with two orbitals per 
donor and acceptor molecule, i. e. two bands per molecular chain, the lower bands being 
partially filled. Applying the recently proposed tight-binding formalism for the dielectric 
matrix in the random phase approximation (RPA) [flif , we determine the hybridized collective 
modes associated with the electron polarization processes, and the corresponding zero-point 
energy contributions to the cohesive energy. We show that the above unusual cohesive and 
optical properties are related through the presence of a soft interband collective mode and its 
hybridization with the intraband charge fluctuations. 

In the tight-binding formalism [Q the matrix elements of the screened Coulomb interaction 
V(r, r') are represented in terms of interband and intraband polarization processes, and the 
corresponding system of RPA Dyson equations reads 

[<W ~ V P*P9 ( q ) n p s ( q ' w )] Vp*Pt ( q ' ^ = V i>'Pf ( q )- W 

Pg 

Here the matrix elements of the bare Coulomb interaction V(v — r') are given by 
Wq)=E e * q(R+e ~ f) I dY I *Vi e (>'-R-e)< / (r'-f)F(r-r')^(r-R-e)^(r'-f). 

R 

(2) 

R are lattice vectors, and e and f are positions of the molecules within the unit cell. The latter 
take values and a/2 for acceptor and donor molecules respectively, with R = N a & + N c c 
for TTF-TCNQ, and R = ^(a/2 + c) + A 2 (a/2 - c) for HMTTF-TCNQ. Here iV a ,c,i,2 are 
integers, and a simplified orthorombic symmetry with one donor and one acceptor chain per 
unit cell is assumed. The wave function tpi e (r — R — e) in eq.^ is Z e -th orbital at the molecular 
site R + e, and p e = (l e ,l' e ) stays for the one-electron transition between l e and l' e orbitals. 
We do not keep negligible contributions due to finite direct overlaps between orbitals from 



neighboring molecules. Having in mind the known electronic spectra of e. g. TCNQ 15 
and TTF |l(| molecules, we specify that the orbitals which form lower (partially filled) and 
upper (empty) bands have the same, A g and i?2u, respective symmetries on both, acceptor 
and donor, sublattices. This means that the interband transitions on both families of chains 
are dipole active, with the dipole matrix elements fi a u\ oriented along the chain direction b. 
With the standard assumption that all products (/?^(r)i/? ; ' (r) are real, we get from cq.[j] the 
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dielectric matrix 



s = 



i-v 0aQa iiQ a -v 0a0d n 0i v 0a i a ii la -v 0aU u ld 

~v 0d0 n 0a i-y 0d0d n 0d -v 0dl n la -v 0dld n ld 

-^i„o a n 0a -^i a o d n 0d i-v la i a n la -v lald n ld 

-v ld0 n 0a -v ld0d u 0d -v ldl n la i-v ldld u u 



(3) 



for the microscopic response to the longitudinal electric field. Here indices a (d) an d l a (<f) 

represent intraband and interband one-electron transitions on the acceptor (donor) chains, 

respectively. IIo a(d) and IIi , d , are intraband and interband bubble polarization diagrams, 

and V's are the corresponding matrix elements of eq. |[ Our main aim is to determine 

collective modes, i. e. zeros of the microscopic dielectric function e(u>) = det[e], in the 

long wavelength limit q — > 0. In this limit the intraband polarization diagrams are given 
■\2, ,2 



by IIo a(d) (q, uj) ~ (qb) uj 2 a ^/ fate b u> , while the interband polarization diagrams reduce 
to rii Q(d) (q, lj) w 2n a ^E a ^/{uj 2 — E 2 ^), after taking into account that the energy gaps 
on both types of chains, E a ^ (« 3eVs) |T^], are considerably higher than the bandwidths 
(~ 0.5eVs) ||. Here = 2 — n a , and uJo a (d) and -E'a(d) are respectively plasmon frequencies 
and differences between upper and lower electron band for the acceptor (donor) sublattices. 
Furthermore, the multipole expansion can be performed for all Coulomb matrix elements 
except for the first neighbor short-range interaction V sr in the a direction. In the q — > limit 
we keep for three types of the Coulomb matrix elements present in eqj^, namely those with 
intraband-intraband, intraband-interband and interband-interband scatterings, the leading, i. 
e. monopole-monopole (Vo e o f — Aire 2 /voq 2 ), monopole-dipole (Vb e i / = iniefif ■ q/vo<Z 2 ) and 
dipole-dipole (Vi e i f = 47r[3(/ix e • q)(^*/ ■ q) — 9 2 A t e ' M/]/3<Z 2 + V sr ), contributions respectively. 
t>o is the volume of the unit cell. Note that the difference between the longitudinal (q||b) and 
transverse (qJ-b) dipole-dipole matrix elements, which will appear in the further analysis, is 
independent of q, as well as of the details of the crystal structure p7t , V\ K \.j, — Vi e \,t — 
(4wii e Hf)/v Q . 

From the form of the matrix (^), one may already recognize the hybridizations which lead 
to collective modes. We start with the upper diagonal 2x2 block which represents intraband 
(metallic) polarizations with long range monopole-monopole interactions. They give rise to 
two plasmon modes, the hybrids of plasmons from acceptor and donor sublattices. For strictly 
one-dimensional electron bands, one of these hybridized modes is acoustic, while the other 
is optic, with a frequency which vanishes in the limit q — > for q_Lb, and is equal to 

The lower diagonal 2x2 block includes interband processes, which, due to finite dipolc 
matrix elements at both types of molecules, induce long-range dipole-dipole interactions. The 
corresponding collective modes 



co 2 



■4±>J (wjj - uo 2 f + 16n a n d E a E d V 2 ald 



(4) 



are hybrids of dipolar modes from two sublattices, tj 2 ^ = E a ^ \E a (d) + ^ n a(d)Vi a{d) i a(d) \ . 
The frequencies of the dipolar modes decrease from u)^i to u^t [i = a, d, ±] as the orientation 
of q changes from q||b to q_Lb p3 [and do not vary as q rotates within the plane (a, c)]. Note 
that the frequency differences Lof^—Lo 2 t are, like the corresponding differences of dipole-dipole 
matrix elements, insensitive to the details of the crystal structure. 

A further hybridization between plasmon u>q and longitudinal dipolar (H) modes takes place 
through the off-diagonal 2x2 blocks in the matrix (ph. Due to the diverging long-range 
monopole-dipole Coulomb matrix elements, the latter do not vanish, provided the lower bands 
are metallic. The frequencies of the three new hybridized longitudinal modes are solutions of 
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the equation 

wV - - u 2 +l ) - uj 2 (uj 2 - uj 2 _ t )(uj 2 - LJ 2 +t ) = 0, (5) 

which we denote by ^1,2, 3- Although the transverse dipolar frequencies u>±t are not affected by 
the monopole-dipole interactions, they enter into eqj^, in agreement with our general results for 
multiband metallic systems |i~4| |. This has an important implication when u>—t, the frequency 
of the lower transverse mode, is small, i. e. when 

y iaU,t ~ — ; =, [p) 

2^n a ndE a E d 

which is, as we argue below, the case in TTF-TCNQ. Then £2i, the frequency of the lowest 
longitudinal mode from eq. |^, is also small and lies below w_t, i. e. 



n 2 



(7) 



The frequencies u) a id)t are presumably comparable with the corresponding values in the ho- 
momolecular crystals of TCNQ and TTF molecules. Dielectric data for TCNQ crystals give 
oj at ~ E a fl9[| . Analogous data for TTF are not available, but, having in mind that the 
spectra of the TTF and TCNQ molecules are similar jl6| , yl| , and that both homomolecular 
crystals have similar cohesion energies of leV |pH , we expect that uidt ~ Ed, too. Thus the 
condition (||) is physically equivalent to the requirement E a ,Ed,Vi a i d ^> Vi a i a ,Vi d i d , which 
in particular means that the dipolar coupling between donor and acceptor sublattices is much 
stronger than that within each sublattice. Since this is not a property of the dipolar long-range 
lattice sums, i. e. the product of the long range parts of V\ a \ a and V\did is of the same order 
of magnitude as the square of the long range part of Vi a id, the dominance of V\ a \d is to be 
sought in the details of the short-range contributions from nearest neighboring pairs of chains. 
We emphasize that, while the soft transverse mode w_t appears due to the strong, mainly 
local, dipolar fluctuations, the accompanying soft longitudinal mode fl\ from eq.^ appears due 
to their coupling of long-range intraband fluctuations. 

The macroscopic dielectric function which corresponds to the matrix (||) is given by p4| ] 

_ (c 2 -Of)(c 2 -^)(^»j) 

£ M — -f—n 2 \f 2 2~\ \° 

for q || b, and em — 1 for qlb. Since these functions correspond to the optical (transverse) 
dielectric functions for the orientation of the electric field parallel and perpendicular to the 
chain direction respectively, we can make a link with experiments, keeping in mind that above 
expressions for cm do not contain contributions from interchain electron hoppings, molecular 
and lattice vibrations (including the 2fc^ CDW collective modes), relaxation processes, etc. 
In particular, we assign the modes Q\ and w_t to the respective excitations at ~ 10 meV 
and ~ 57 meV, observed at 100 K in the infra-red measurements |ITJ, [TT|, and the mode SI2 
to the excitation at ~ 0.75 eV, originally interpreted as a Drude edge~|| [TJ (fig. [j]). The 
modes u) +t and ^3 are most likely situated in the frequency range of few eVs, not investigated 
experimentally in detail and usually attributed to the range of interband transitions. It is 
noteworthy to mention in this respect a strong anomaly at ~ 4.5 eV observed in the early 
reflectivity data M. 

The two most interesting modes, denoted here by Oi and w_t, invited various theoretical 
explanations. Here we list arguments in favor of the present one. First, our model predicts 
the optical activity in the chain b direction, in agreement with experiments |T7J, [HJ, and in 
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Fig. 1. - The real part of the dielectric function for TTF-TCNQ at 100 K with the data taken from 
Refs. |ll|] (triangles) and [hoi (circles). The full line is Re em (eq. fef) with the zeros at fh = 10 meV, 
= (30 meV, and SI3 = 6eV, and poles at ui- t = (57 + i42)mev and £j+t=4.5 eV. The imaginary 
part in ui-t is chosen to fit the experimental width of the pole in Re em. 



contrast to the interpretation in terms of lattice vibrations |Tg . The latter are optically active 
only when polarized along the transverse a direction, in which the oppositely charged chains 
are aligned. Our suggestion that the regime (||) is appropriate for TTF-TCNQ is based on the 
observation that the frequencies Oi, uj-t, and SI2 are much lower than those of other modes 
from eq.^|, which lie in the range of bare interband dipolar modes (~ 3 eV). Furthermore, 
the experimentally determined difference of ~ 0.7 eV between frequencies H2 and u~t ju|, 
which are the longitudinal and transverse edges of one hybridized branch (out of three), is 
comparable to the difference of ~ 1.5 eV between analogous frequencies in, e.g., the TCNQ 
crystal Q (although the characteristic values of latter are of the order of a few eV). Since 
these differences measure the differences between longitudinal and transverse lattice sums in 
the matrix elements Vn, it follows that the long-range contribution in Vi a % d is comparable 
to those in Vi a i a and Vi d i d , in accordance with the already mentioned general result Q. In 
other words, the small values of both cu_ t and VL 2 may originate only from the short-range 
q-independent part of V\ a \ d , in agreement with our arguments below eq.^. 

The onset of CDW order below ~ 50K , which causes the appearance of the well-known 
features in the far infra-red range due to the CDW dynamics jjJJ , have also important effects 
on the modes of eq|| At first, the low- lying longitudinal mode Oi is eliminated |lj, [f| due 
to the formation of the insulating gap Al at the Fermi level, in agreement with our general 
predictions JlJ]. The modulation of the crystal structure in the insulating phase below the 
temperature of the Peierls phase transition also causes a weak splitting of remaining transverse 
and longitudinal modes , in analogy with Davydov splitting of molecular excitons . Such 
splitting is indeed observed for the mode w_ t at ~ 50 meV fll3[| . We expect that more refined 
measurements would show the same splitting for other modes, e.g for the mode SI2 at 0.75eVs. 

It remains to consider the contribution of the polarization modes (|5|) to the cohesion. 
Extending the standard RPA procedure to the present multi-band system, one arrives at 
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the cohesion energy per unit cell 

E coh = Ein a - E M n 2 a + ^ ~ E ° ~ Ed > ^ 

where uij are hybridized modes from eqsjiJJ^. Due to the condition^]), the expression (|^) 
reduces for n a < 1 to 

E coh » -(2C -D- Ej)n a + (C - ~ - ~ - E M )n\ - D, (10) 

Here A(D) = £ q K w i* w / 2NE M> C = £ q V? aU /N(E a + E d ), and E 1 and E M are 
afhnity-ionicity and Madelung energy respectively Noting that, in accordance with foregoing 
arguments, A, D <C C, it follows from eq.[l(] that only with 2C — Ej > and C — Em > is 
partial charge transfer possible. In the case of TTF-TCNQ this means that, roughly, C > 2eV, 
which just compensates the "Madelung defect" on one side, and is consistent with the condition 
for the u)—t mode softening ([|) on the other. As was already pointed out, the dominant 
contribution to C is expected to come from the nearest-neighbor interaction V sr between 
acceptor and donor pairs, i. e. from the a direction which is indeed the direction of smallest 
compressibility in TTF-TCNQ Q. This interchain coupling is expected to be responsible for 
the fractional charge transfer in other OCCs as well. The fractional charge transfer is thus 
predicted to be relatively insensitive to the type of stacking in the c direction, which explains 
the similarity of the observed charge transfers in e. g. TTF-TCNQ and HMTTF-TCNQ. 

Let us now turn our attention to the relative stability of TTF-TCNQ and HMTTF-TCNQ 
lattices. Note that the short-range contribution from nearest neighbors in the a direction to 
the sum C does not enter into the difference of corresponding cohesion energies (jiol), while 
the contributions to the sums A, D, C from the neighboring molecular pairs in the c-direction 
which now become relevant, can be represented in the Friedel dipole-dipole form. Taking for 
simplicity n a = 1, we come to the condition that TTF-TCNQ lattice is more stable than the 
HMTTF-TCNQ lattice provided that 



E a VElE2{E a + E d ) 



1 1 



6 [(a/2) 2 4-^13 



>E M (H)-E M (T). 



(11) 

E M (H) and E M (T) are Madelung energies of HMTTF-TCNQ and TTF-TCNQ lattices re- 
spectively. The left-hand side of this expression amounts to the second order van der Waals 
contribution. It is a slight generalization of the original Friedel criterion, since E a ^ E^ 
and the contribution from the diagonal a/2 + c direction is taken into account. We note 
that the characteristic energies which are responsible for the relative stability of TTF-TCNQ 
and HMTTF-TCNQ lattices are much smaller than those which determine the total cohesion, 
including the contributions from the coupling in the a direction which is decisive for the 
fractional charge transfer. Indeed, the simple calculations of the differences of the Madelung 
and the van der Walls energies in the point charge approximation by using the numerical 
values for the molecular polarizabilities , suggest that both differences are of the order 

of O.leV per unit cell. 

In summary, we have shown that the interband collective modes are responsible for the 
crystal stability of OCCs, while the hybridization of of these modes with intraband plasmons 
explains the low frequency optical data in TTF-TCNQ in particular. We argue, and this is 
the crucial point of our model, that the frequencies of three collective modes are rather low, 
i.e. below 1 eV. It is worthy to note that, since one of them falls into the frequency range 
of the CDW instabilities (< 10 meV), we are confronted with the interesting possibility of a 
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mixing of these two types of correlation at low temperatures. This question is left for future 
investigations. Clearly, the concepts developed and successfully used here can also be applied 
to the other metallic systems with large molecular polarizabilities. 
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